Effective potentials for 6-coordinated Boron: a structural approach 
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We have built an ab-initio LDA energy database with over 
60 hypothetical extended structures of pure Boron, in each of 
which the coordination environment of each atom is equiv- 
alent. Focusing on eleven 6-coordinated structures, which 
are most relevant to the observed Boron phases, we deduce 
a nearest-neighbor-only 2-body interaction, refined by poten- 
tials that depend on angles and asymmetry within the local 
atomic environment. The resulting effective potentials de- 
scribe the 6-coordinated structures with an average error of 
0.1 eV/atom, and favor the environment most often seen in 
real Boron. 

PACS. 34.20.Cf, 61.50.Lt, 71.20.Mq, 33.15.Dj 



Boron as the prime candidate for a covalent quasicrys- 
tal has inspired activities in both experimental search 
and theoretical assessment, leading to identification of a 
new phase |l|] and interests in nanotubes ||. Accurate 
and fast structural energy estimators, capable of handling 
complex extended structures, can facilitate further inves- 
tigations of these proposed structures, (through Monte- 
Carlo or Molecular Dynamics exploration), and might 
elucidate the structures of the unsolved or amorphous 
phases. While ab-initio methods have been used for com- 
paring specific Boron structures , we seek a less com- 
putationally intensive effective theory, in the style of the 
environment-dependent effective interatomic potentials, 
which have successfully described Silicon bonding 
but are not yet available for Boron. 

The bonding behavior of Boron is unique among the 
elements and allows a number of complex structures. 
The atomic structures have been refined for only the 
metastable au,Tc,o and the stable /3ios @, while various 
complex phases remain unsolved ||. Icosahedral clus- 
ters of various scales are present in all the known atomic 
arrangements, leading to unusual coordination environ- 
ments and geometric patterns Q, quite unlike those in 
the well-studied covalent (e.g. Silicon) systems. In the 
case of the structurally simplest a±2 phase, accurate elec- 
tronic band structure ||, lattice dynamics studies and 
force models jlO| all confirm the presence of 3-centered 
bonds, and the coexistence of soft metallic and strong 
covalent bonds both within and between icosahedra [^l) . 

To design a classical potential (total energy function) 
that describes not only the uncommon bonding mech- 
anisms and large-scale structural features, but also the 
unforeseen local configurations that may arise in qua- 
sicrystal models, we assume a local form E to tai = E i> 
where Ei, the site energy of each atom, depends on its 



local atomic environment, as parametrized by bond dis- 
tance, angles, and quantities involving more atoms. 

Featuring the coordination number Z in environment- 
dependent energy functionals for Silicon, Tersoff J^] pro- 
poses ^-dependent bonds, while Bazant et al || formu- 
late ^-dependent 2- and 3-body potentials to better cap- 
ture the angular behavior. Making no assumptions for 
Boron, we investigate this type of description in its most 
general form of a independent N-body expansion, 
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where j,k ^ i labels neighbors to atom i, r the magnitude 
of the interatomic distance, and 6^ the angle formed by 
two neighbors j, k and center i. 

In this letter we discuss a structural approach that 
naturally determines such potentials and allows straight 
forward assessment of the angular forces without assum- 
ing an analytic form. With an extensive energy database 
of the so-called "uniform" structures, constructed to have 
nearly identical bond lengths but large angular varia- 
tion, we sample bonding mechanisms in widely different 
topologies and configurations. Focusing on 6-coordinated 
(Z6) structures, which are relevant to the real phases, we 
determine that their cohesive energies can be sufficiently 
described by Z6 potentials of the form 

E ™ = £/(^)+£Kwi^§)+M&) ( 2 ) 

3 3+k 

where j, k are now restricted to "nearest" neighbors of i, 
and the last term is a function of the site asymmetry £j 
to be defined later. 

In a "uniform" structure Jl2], all sites are geometri- 
cally equivalent by space-group symmetry operations, so 
that there is only one type of atomic environment oc- 
curing throughout the entire structure. We investigate 
those uniform structures in which all the nearest neigh- 
bor distances are equal (or nearly equal), so that angular 
effects are more clearly detected. Our database of over 60 
uniform structures Jk|, spanning coordination numbers 
Z from 3 to 12, and exhibiting large angular variations 
for each Z, serves as a set of simple but topologically dis- 
tinct representatives for different regions of configuration 
space. In fact, Boron compounds adopt certain of these 
structures as the sub- lattice for the Boron atoms p3|. 

ZQ coordinations are particularly important in Boron 
bonding, as they show up in 80% of the sites in /?io5, 
and 50% in a\2- Therefore we focus now on the uni- 
form Z6 structures, their bonding, and the effective po- 
tentials determined from them. Other coordinations are 
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discussed in In our list of Z6 structures (see Fig- 

ure ||), besides the simple cubic (SC) and the triangular 
(tri) lattices, there are two 3D tilings: the simple cu- 
bic array of cuboctahedra (SC-co) and the diamond ar- 
ray of vertex-joint tetrahcdra (T-lattice) From two 
Z4 planar lattices, the 3.4.6.4 ("i-dodec," for resembling 
a set of inter-penetrating dodecagons) and the 3.6.3.6 
(Kagome), we get Z6 structures by stacking them with 
a separation of unit length (defined as the nearest neigh- 
bor distance.) There is also the Primitive Orthorhombic 
packing of icosahedra (PO-ico) [|l6| with inter- icosahedral 
distance set at unity. 

We have constructed some novel Z6 structures as well. 
When we pucker a planar lattice, and mirror the opera- 
tion in the next layer, each vertex gains one additional 
neighbor from either the layer above or below. So from 
the Z5 lattices 3 2 .4.3.4 (called a) and the 3 3 .4 2 (called 
//), we get the Pa and P/i Z6 structures, ("P" for puck- 
ered). When we roll the triangular lattice into a single 
infinite tube with a circumference of 8 edges, we get the 
"tube-tri" structure [§). Again, by rolling the 3 2 .4.3.4 
(cr) lattice into a tube with a circumference shown by 
the dotted line in Figure [l]-i, and arranging such tubes 
symmetrically in a square array (view along tube axes in 
Figure |l|-j), we get the "tube-cr" structure, in which each 
atom has 5 neighbors on its own tube, and one more from 
a neighboring tube. 

We show in Table | the distribution of "geometric 
quantities" found in the coordination environments of 
our Z6 structures, which will become helpful informa- 
tion in potential fitting. These consist of N 2 and N 3 , the 
number of neighbors on the second and third shell, (at 
distances of y/2 and y/H, respectively), and frequencies 
of various "two-neighbor" angles, (which are formed at 
a central site and point toward two nearest neighbors). 
We also describe the asymmetry of these atomic environ- 
ments, parametrized by 



3 V Jl 



(3) 



where j labels the nearest neighbors, and Vj the vector 
that points from the central atomic site to j. The denom- 
inator is just the average distance to nearest neighbors. 

We perform ab-initio total energy calculation for each 
structure in our database |17|, in the local density ap- 
proximation [ ]l8{ with extended norm and hardness con- 
serving pseudo-potentials JllJ. 

For each Z6 structure, the edge length "i?" represent- 
ing the nearest neighbor distance may be varied. In Fig- 
ure ||, we present the cohesive energy (per atom) of the 
structure as a function of R in a range of physical inter- 
est. In the inset, where R scales from 1 — 5A, we show 
the resulting single-well curves for three representative 
structures, SC, tri, and PO-ico PC|| . The energy curves 
of i-dodec and Kagome fall closely onto the SC curve, 



and similarly the T-lattice onto the tri, and SC-co onto 
the PO-ico. Beyond R = 3A, all curves converge, and 
saturate toward the non-bonding limit. 

In the main figure, we enlarge the [1.6, 2.0] A region, 
(within which bond lengths from different structures all 
fall ) , to clarify the difference among the wells of the 
energy curves. Our lowest energy Z6 structure is the PO- 
ico, a topologically different icosahedral packing. Within 
the next low energy structure, SC-co, is the cuboctahe- 
dron that has been compared to a distorted icosahedron 
using molecular orbitals theory pl[ |. For comparison, we 
include the lowest energy curve of the "ideal-ai2" phase, 
topologically same as a± 2 but with all nearest neighbor 
distances set equal. 

To construct an effective theory, it is natural to start 
with a typical single-well 2-body potential, with nearest 
neighbors close to the minimum of the well, and further 
neighbors contributing at the tail of the well. This will 
not work for our Z6 energy curves. The large variation 
of N 2 and N 3 in Table | implies that the energy contri- 
bution from the further neighbors will lead to a large en- 
ergy variation, on the order of the well-depth, among the 
structures. This contradicts the small energy differences 
among our energy curves in Figure |^. Indeed, numerical 
fittings with several single- well forms all substantiate this 
argument. We conclude that the 2-body potential has to 
be independent of N 2 and N$, i.e. it must be essentially 
truncated after nearest-neighbor range. 

Since we have shown the important result that the first 
term in Eq. |l| sums only the nearest neighbor j, our Z6 
structures must have basically the same 2-body contribu- 
tion of 6/(i?) |Q. The energy differences among these 
structures would have to come from higher-order poten- 
tials. We choose the simple cubic energy curve to serve as 
the arbitrary reference in taking energy differences, and 
to give the function f(x). We propose a 3-body poten- 
tial that depends only on the two-neighbor angles. We 
also include a term h(£) (compare Eq. ||) linear in £ to 
test the degree of energy dependence on the asymmetry, 
which involves all Z neighbors. Allowing for dependence 
on the scale R, the cohesive energies (per atom) would 
lead to a set of linear relationships, 



AE rn (R) = )J N m {9 a )g(R, R, 9 a ) + c(R)^ r , 
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where A.E m is the energy difference between the m-th 
structure and the reference, a labels the two-neighbor 
angles, N m the number function per atom, and c the 
linear coefficient of £ that depends on R. 

A major problem in determining the functional form 
of the angular potential g(R,R,0) is that the frequency 
of a particular angle 6 a cannot be varied independently 
of 6 a i to display its individual effect on the cohesive en- 
ergy. Normally this problem is avoided by fitting to an 
assumed form j(| , but for Boron, where angular behavior 
is still being explored, we have no analytic guidance as 
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to the form of this potential. However, the small number 
of distinct angles 9 a appearing among the structures in 
Table | helps to resolve this. 

Since all nearest neighbor lengths are nearly equal, we 
can separate the independence from the angular effects. 
The low number of distinct angles that appear in our 
Z6 structures allows us to determine the shapes of the 
angular potential without further assumptions. With 
the eleven Z6 structures and the data from Table |, 
we can solve numerically, in the least-squares sense, for 
g(R, R,9 a ) and c(R), not assuming any functional form. 
Doing this at different R gives the radial dependence of g 
and c. The resulting fit, shown in Figure 0, suggests that 
the angular potential is linearly monotonically increasing. 
The independence for each angle and c(R) shows rapid 
decay. Similar potentials are found for Z5 structures jl3) . 

Although the angular potential is obtained from the 
correlation of the number of angles with total energy, and 
not from measuring the direct effects each angle has on 
the site energy, this result is by no means trivial: for the 
potential value at each of the 7 angles is an independent 
parameter in the fit, yet the resulting g(R,R,6) looks 
smooth as a function of (R,9). Furthermore, we find 
that a separable form, i.e. g{rx,r 2 ,9) = g(ri)g(r 2 )A(9), 
is sufficiently determined, and a good approximation. 

When we apply our potentials g and c to the Z6 
"inverted-umbrella" environment |lC| ] as found in the 
ideal-ai2 structure, its site energy is 1.0 eV lower than 
that of the second best environment (PO-ico) among our 
structures. In fact, a Monte Carlo search among ZQ en- 
vironments having all two-neighbor angles greater than 
60°found no lower energy configurations. Since the effec- 
tive potentials were constructed from environments all 
higher in energy, it is encouraging that they predict the 
realistic "inverted-umbrella" to be lowest in energy. Fur- 
thermore, our potential disfavoring 180° is in agreement 
with the tendency for Boron clusters to buckle in a hexag- 
onally coordinated environment 

The angular and asymmetry potentials constructed so 
far can account for most of the 1-2 eV/atom energy vari- 
ations among the ZQ structures, with an average error of 
0.14 eV/atom in the bonding range. The energy ordering 
of our Z6 structures are mostly retained by the poten- 
tials. Although we still need to treat other Z's, and test 
cases when not all neighbors are at the same distance, 
we see a promising start in the local potential descrip- 
tion for the general Boron system, without many-body 
terms involving, e.g. all the atoms in an icosahedron. 

We thank M. Teter, D. Allen, J. Charlesworth for 
providing code and support in the LDA calculation, A. 
Quandt, M. Sadd for comments, K. Shirai, P. Kroll and 
R. Hoffman for discussion. This work is supported by 
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TABLE I. Local environments of Z6 structures. 



Structure N2 


N 3 


60° 


90° 


108° 


N{0 a ) 
120° 135° 


150° 


180° 




SC 


12 


8 





12 














3 





tri 





6 


6 








6 








3 





SC-co 


7 


6 


2 


7 





2 


4 








0.586 


T-lattice 





12 


6 








6 








3 





i-dodec 


10 


6 


1 


10 





1 





2 


1 


0.732 


Kagome 


8 


4 


2 


8 





2 








3 





PO-ico 


1 


10 a 


5 


1 


7 








2 b 





1.854 


Pa c 


4 b 


12 a 


3 


3 b 


4 b 


l b 


2 b 


2 b 





0.263 


P M 


5 b 


10 a 


3 


4 


3 b 


2 


2 b 





1 


0.771 


tube-tri 





6 a 


6 








6 b 


1 





2 d 


1.163 


tube-cr 


5 b 


It3 


3 


5 b 


2 b 


2 





3 b 





0.936 



a Neighbors between the second and third shells, occuring at 
various distances, are binned in A3. 

b These bins have tolerance of 0.14 unit distance, or 6°in angle. 
c Nearest neighbors not all at unit distance pq ], 
d This is at 169°. 
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I ■ I ( j . 1. The less familiar Z6 structures, with Bravais lattice 
indicated by *. For structures in layer form, • are nearest 
to the viewer, and o furthest, (a) SC-co (• • o stacking), 
(b) T-lattice, in which the Kagome layer has an (ABC) type 
stacking, (c) i-dodec, (d) Kagome, (e) PO-ico, (f) Pa, (g) Pfj,, 
(h) tube-tri, (i)&(j) tube-cr (see text). 
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FIG. 2. Cohesive energy Ei(R) for Z6 structures. In the 
inset, the three representative structures are SC (line), tri 
(short-dash), and PO-ico (long-dash). The region showing the 
energy wells is magnified in the main figure. The structures 
are, in order of lowest to highest energy at R — 1.8A, ideal-ai2 
(line), PO-ico (long-dash), SC-co (□), tube-cr (o), T-lattice 
(line), P-/i (long-dash), tri (short-dash), P-a (x), tube-tri 
(line), SC (line), i-dodec (o), and Kagome (x). 
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FIG. 3. Angular and asymmetry potentials. The angles 
take on values from 60°to 180°. Data in the "Asym" column 
signifies the coefficient c(R) in the asymmetry potential. The 
g(R,R,6) and c(R) are illustrated by symbols {*, o, x,+,o}, 
corresponding to R = {1.6, 1.7, 1.8, 1.9, 2.0}A. 
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